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1 Introduction 



The investigation of two-dimensional classical and quantum integrable field theories with one or 
two boundaries is now several decades old (see for example [I]- [9]), yet it is incomplete with many 
interesting outstanding questions remaining to be answered, even in the context of the much- 
studied affine Toda models. More recently, there has been some interest in the definition and 
properties of integrable defects (see, for example [ID]- [22]), which is really the study of allowable 
discontinuities, in the sense of permitting field discontinuities that do not upset integrability, 
and already there have been suggestions [T5| [T8| [20] concerning possible relationships between 
the two ideas. The purpose of this article is to add to these suggestions having in mind models 
that allow solitons and using the sine-Gordon model as the prime example. 

Basically, a defect is a bulk phenomenon but placing a defect close to, or even at, a boundary 
will effectively modify the reflection matrix associated with the boundary in several ways. First, 
the defect is able to introduce additional parameters, and second, even if in the absence of the 
defect a reflection matrix contained no explicit dependence on the topological charge residing 
at the boundary - since boundary conditions only exceptionally preserve topological charges 
carried by solitons - the defect will (almost) inevitably introduce dependence on topological 
charge in the modified reflection matrix. So, it would seem sensible to reconsider the boundary 
Yang-Baxter equation, or reflection equation, allowing the boundary to carry topological charge, 
and the reflection matrix to be infinite-dimensional, then attempt to classify its solutions. 

If the sine-Gordon model is considered as the primary example, on the grounds that it is the 
simplest of the affine Toda field theories, one might expect the solutions to the refiection equa- 
tion to be one of three types: (1) the coefficients of the reflection matrix could be independent 
of topological charge (basically, this is the Ghoshal-Zamolodchikov solution though reformu- 
lated slightly); (2) the reflection matrix is related to a Ghoshal-Zamolodchikov solution modifled 
by a defect; or (3) neither of these - in which case the reflection matrix may describe a new 
type of boundary condition, which might turn out to be a defect 'fused' with a boundary. The 
latter idea might appear strange at flrst sight yet it is worth recalling that defects themselves 
may be fused together and the essential difference between an unfused and a fused pair is the 
evident fact that a soliton can propagate in the gap between the unfused pair and there is no 
gap between a fused pair |22]. One need look no further than a Neumann boundary condition 
(placing the boundary at a; = 0), Ux{0, t) = 0, to realise that a boundary can store an unlimited 
quantity of topological charge since N solitons approaching the boundary eventually reflect as 
N anti-solitons, or vice versa, with the boundary picking up 2N units of charge. In effect, the 
boundary value of u changes by AttN during the process. On the other hand, a Dirichlet bound- 
ary condition, for example u{0,t) = uq, ensures solitons reflect as solitons and the boundary 
charge does not change during the process. 

Under reasonably general assumptions, effectively meaning the absence of any dynamical vari- 
ables conflned to the boundary, classically integrable boundary conditions have been determined 
[5] for any of the (real) affine Toda field theories in the following form 



where m is a multicomponent field, the vectors are the Euclidean parts of the affine roots 




(1.1) 
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described by any of the Dynkin-Kac diagrams [26], the coefficients br are constants. These 
boundary conditions possess curious properties [6]. With a single exception, the most unexpected 
concerns the affine Toda field theories based on an\ dn\ root data. Apart from the 
sine-Gordon model (the exception, based on a^^^ data), which has two free parameters 61,62 
associated with a general integrable boundary condition, these models do not appear to allow 
any free parameters associated with a boundary. Instead, there is a discrete set of non-zero 
values for the coefficients br amounting to a set of sign choices (Neumann is always a possibility 
with br = for all r). The other models (based on non-simply-laced data) fare a little better, 
in the sense that there are usually one or two free parameters but never as many as the rank 
of the associated affine algebra. In the complex affine Toda models (where there are complex 
solitons), there are other boundary conditions [7] of the form 

lmu(0,t) = 27rA, Reu^(0,t) = 0, (1.2) 

where A is a weight. Clearly these are a mixture of Dirichlet and Neumann conditions, and they 
are soliton-preserving. In the complex models the conditions fll.ip are not soliton-preserving 
and will inevitably change the topological charge on the boundary. It is possible that the 
assumptions made to reach these results are too stringent but it has not yet been found how 
to relax them. It is also possible that dressing a boundary with one or more defects may 
add additional parameters and provide a clue towards the goal of generalising the boundary 
conditions. On the other hand, this is easier said than done because, so far, defects are only 
established for the affine Toda models based on an\ n > 1 and af^ [211 [23]. 



2 Generalised reflection matrices 



Consider the sine-Gordon model, a typical refiection matrix is 

KiiO), (2.1) 

where the labels a, 6 = ±1 (or simply ±) represent the incoming and outgoing soliton, the 
labels a, 13 represent the initial and final topological charge carried by the boundary, and 9 is 
the rapidity of the incoming soliton. In much of the earlier literature the boundary topological 
charge labels are simply suppressed but here it will be useful to keep them explicit. Clearly, 
topological charge conservation requires 

a + a = b + P, (2.2) 

and the boundary labels must be either even or odd since they may only change by ±2. The 
refiection matrix will also depend upon the bulk coupling and boundary parameters. In the more 
general affine Toda models, both sets of labels should be replaced by suitable sets of weights. 

The refiection equation, or boundary Yang-Baxter equation [1], is a condition that ensures 
compatibility with the bulk soliton S-matrix. Thus 

RfM s^e, + e^)R;iie,) si:{e, - e,) = sii{e, - e,)Ri%{e,)si:{e, + e,)Rl}{e,), (2.3) 
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where repeated indices are summed. In the sine-Gordon model, the bulk S-matrix [27] is given 
by a 4 X 4 matrix depending on the parameter G = G12 = {61 — 62) whose elements different 
from zero are: 
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p,(e) = a(e)p,(e), 
p,(e) = 6(e)p,(e), 



(2.4) 



with 
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The multiplicative factor Ps{Q) is given by: 



p,(e) = I<l^i^2Ml±i)nfi.(e)fl.(«-e). . = ^5. (2,5) 

with 



/c=l 



^ r(z + 2A;7)r(l + z + 2A;7) 

^ r(z + (2fc + l)7)r(l + z + (2A;-l)7)- 

There are many solutions to (12. 3p that may be written 

where the reflection factor is written in block form, and it is understood the blocks are labeled 
by the solitons. Typically, the coefficients appearing in (12. 6p will be functions of a, the initial 
charge on the boundary, besides the rapidity, bulk coupling and boundary parameters. However, 
there is a particular two-parameter class of solutions, referred to as the Ghoshal-Zamolodchikov 
solution, where there is no dependence on a other than via the Kronecker deltas. In this case the 
latter are there merely to keep track of the topological charge. This solution has the following 
form, 

where ri, r2, /cq, ^0 are constants and x = e'^ . Without losing generality one may impose rir2 = 
1, ko = lo, hj removing an overall constant factor and using a similarity transformation. The 
requirements of crossing and unitarity restrict the overall scalar factor a{d) [1]. 

A more general solution to (12.31) can be obtained quite easily by assuming the elements s±(a, x) 
remain proportional to — and the diagonal elements are proportional to a cubic in x and 
1/x. Explicitly, apart from an overall function of rapidity, the solution is: 

r+(a;,x) = (x^ — l/x^) (^r^q'^'^^x — r4g~°~^/x) -|- rix + r2/x, 
r_(a, x) = (x^ — (r4g~""'"^x — r-sq'^"^ /x^ + r2X + ri/x, 

s+(a, x) = {x^ - 1/x^) {ko + + fcag""), s_(a, x) = {x^ - 1/x^) (/q + /ig" + hq"") 
= -rl, k2l2 = -rj, kilo + q^koh = grsrg, kok + g^/cs/o = qrir^. (2.8) 
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Clearly, the expression (12 .Tp lies within this class of solution and is obtained by setting = 
= = ki = k2 = h = h- Another interesting special case concerns the choice ri = r2 = for 
which the quantity (x^ — becomes an overall factor. Actually, assuming all terms in (12. 6p 

are polynomial in x, it is not difficult to show that the polynomials must have common zeros, 
except for the roots of = 1, which may be present in the off-diagonal terms without also 
being roots of the diagonal terms. For this reason, the solution given above is general. More 
detail concerning these solutions is supplied in appendix A. 

If a defect is placed before the boundary this will modify a reflection factor though not merely 
by altering its dependence on rapidity. Since a defect also carries topological charge a new set 
of labels will be added to the modified refiection matrix H^]. Thus, in detail, 

Rliio) = T:im'i{o)f',^{oi (2.9) 

where a,/?, 7 are topological charges associated with the defect, T represents the transmission 
matrix for the soliton approaching the boundary, while T represents the transmission matrix 
for the soliton after refiection from the boundary. In fact T{9) is the inverse of T{—9). The 
combination on the right hand side of (12. 9p automatically satisfies a suitable generalisation of 
the compatibility relations (12.31) . provided the transmission matrix is also compatible with the 
bulk soliton S-matrix. The latter requires, 

ST^{ea - 9,) TZie^) T2^{e,) = rfje,) T^j{e^) s:^,S0a - e,). (2.10) 

If the refiection matrix on the right hand side of (12.90 is proportional to 5^ then the new refiection 
matrix on the left hand side is also proportional to 5f and satisfies precisely (12. 3p . 

In the sine-Gordon model solutions to (I2.10p are known and may therefore be used to construct 
solutions to (12. 3p . Owing to topological charge conservation, they all possess the general shape 

TbP(Q\ - ( ^^^^^^^^^ K«,a;)5f-2\ 

T„J^)-p(^)^^^^^^)^,+2 )■ (2-11) 

In [22], the Konik-LeClair solution [11] to (I2.10p was generalised to 

r^P(n^ _ ( («+^^"^'^^' + «-9"/'a;) 5f /.(«) 5^^ \ 

with 

/i(a) A(a + 2) — /i(a — 2) A(a) = (g — q~^) {a_d-q°' — a+d^q~°'), 

which implies 

/i(a - 2) A(a) = a.d.q''-^ + a+d+q-''^^ + 7. (2.13) 

Actually, it has not proved possible to generalise (12.121) further and it is most probably already 
the general solution (see also [28]). The constraint (I2.13P is satisfied by choosing, for instance, 

/i(a) = 6+g-°/2 + 5_g"/2^ A(a) = c+g-°/2 + c_g°/^ a± rf± - 6± c± = 0, 

where a±, 6±, c±, d± are otherwise free (complex) constants. Then, setting Q = q~°'^'^, the 
expression (I2.12p becomes 



(2.14) 
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Its inverse is given by: 
where 

A = a^d^qx^ + a_,.(i_g~"'^x~^ — hj^c^q^^ — b_c^q. 

Take the reflection matrix appearing on the right hand side of (12. 9p to be diagonal, for example 
a Ghoshal-Zamolodchikov solution (12. 7p with ri = l/r2 = r, fco = /o = 0, corresponding to a 
Dirichlet boundary condition, 

Then, the result of the calculation fl2.9|] has the general shape fl2.6p with coefficients of the 
following form (up to an overall factor) 

r+(a,x) = (x^ — {a^d^q"'~^^rx — a^d+q~°'~^r~^x~^) + x {gr — hr~^) + x^^ {gr~^ — hr) , 
(x^ — {a+d^q~'^~^^r~^x — a-d^q°'~^rx~^) + x {gr^^ — hr) + x~^ {gr — hr~^) , 



s+(a,x) = (x^ — X ^) (a+6_r -"^ — a^b^r — a^h^q^r + a+b+q 



s_(a, x) = (x^ — X ^) ((i_c+r — (i+c_r -"^ + (i_c_g°r — d+c+g °r -'^) , (2-17) 

with 

g = a_d+q~^ + a+d_q, h = b^c^q"^ + 6+ c_g. 

The Konik-LeClair [11] (type I) solutions are a particular case of (I2.14p . Set, for example, 
a_ = (i+ = c_ = 6+ = 0, multiply (I2.14p by x and collect an overall factor. Then, (12.141) 
becomes 



with u^^"^ = (d^/a^Y^"^ , b^ = c+, e = b-(a^d^)~^^^ and pj is an overall factor constrained 
by the usual requirements of crossing and unitarity. By performing a similarity transformation 
the a dependence of the off diagonal entries can be eliminated and the type I defect matrix, as 
written in [22], is recovered. Then, the corresponding coefficients of the R matrix f)2.9p are: 



r_|_ = X (r g — e^r^^q^^) + x^^ (r^^q — e^rq"^) , 
r_ = X {r~^q — e^rq"^) + x^^ (rg — e'^r~^q~^) , 
s_i_ = (x^ — x"^) er'^v'^/"^, 



s_ = (a;2 - 5ri/i/^ (2.18) 
which is equivalent to the full Ghoshal-Zamolodchikov solution. 

On the other hand, a 'type IF defect is obtained from (I2.14p by setting a+ = = 1, a_ = 
—6+ (i+ = —6+ b- g^, 0+ = — 6_ g^, c_ = —6+ and multiplying by x. Then the expression in 
[22] is recovered. That is, 

T ^/^rm-. rm^ (Q" - 6_ 6+ x^) 5f x {b^Q- + b.Q-^) 5^' \ 



where the form of the scalar function p// is also constrained by crossing and unitarity. The 
corresponding coefficients of the reflection matrix (12. 9p are: 

r+ = (x^ — — 6+ 6_(5~^"~^rx) + x {gr — hr~^) + {gr~^ — hr) , 

r_ = (x^ — (— + 6+6_(5~^""''^rx~^) + x (fir"^ — hr) + x"-*^ ((yT — hr~^) , 
s+ = (x^ - x"2) (6_r-i + \h+\%^r + h+Q^°'r~^ + bH+Q~'^°'r) , 

S„ = (x2 - X-2) (6„g-V + |6+|26_g-V-l + 6n+g2a-8^-l ^ 5^g-2a^^ ^ (2.19) 

with 

g = (|6+n6-P + 1) g-^ = - {\b^\'Q-' + g-^. 

These entries correspond to a new solution of the reflection equation (12. 3p . in which the explicit 
dependence on the label a cannot be removed by a similarity transformation. 

It is tempting to suggest, given the Lagrangian description of a type II defect [2T1 that the 
system associated with this reflection matrix is described classically (dropping explicit reference 
to mass scale and the bulk coupling constant), by the following Lagrangian density 

£b(m. A) = e{-x) Cu + S{x){u\t - B{u, A)), = u^j2 - ul/2 - (e" + e"") , (2.20) 

where u is the sine/ sinh- Gordon fleld (depending on whether u is imaginary or real), A is a time- 
dependent fleld deflned at the boundary and B{u, A) is a functional to be determined. Forcing 
the energy-like spin three charge to be conserved (using the argument appearing flrst in [1]), the 
form of the functional B is constrained to be: 

B{u,\)=e'/'f{u) + e~>^/'g{u), (2.21) 

with 

f{u)g{u) = /i+e"/2 + h^e-""/^ + 2 (e" + u"") + ho, (2.22) 

where h±, h^ are free constant parameters. Since redeflning A — A+A(m) changes the Lagrangian 
by a total derivative, it follows that the boundary part of the Lagrangian density is essentially 
a three parameter expression, which can be represented conveniently by (12.201) and (I2.2ip on 
choosing 

f{u) = /o + y2(6e"/2 + rie"/2), g[u) = go + V2{b-^e''/^ + he""'^). (2.23) 

Notice that the most general solution of (I2.22p has an additional parameter. However, when 
the expressions of the functionals / and g are used in (I2.2ip one parameter can be eliminated 
by rescaling the fleld A. It can be verifled that this expression for the boundary potential 
corresponds to the expression obtained by adding together a type II defect potential with the 
Ghoshal-Zamolodchikov Dirichlet boundary potential (see [18], [22] for similar arguments). 



3 Discussion 

Some comments are in order at this point. A defect can be used to create solutions to the 
reflection Yang-Baxter equation (12. 3p via (12. 9p . Using the reflection matrix (12.160 representing 
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any Dirichlet boundary condition and the transmission matrix of a type I defect the Ghoshal- 
Zamolodchikov solution (12. 7p is recovered in f l2.18p . In this case the reflection matrix has no 
dependence on topological charge other than via the Kronecker deltas. Note, the simplest 
Dirichlet boundary condition imposes u{0,t) = 0, and its reflection factor has the form fl2.16p 
with r = 1; it is proportional to the identity. This particular boundary condition preserves the 
bulk symmetry u — )■ —u and, using a type I defect, the solution of (12. 3p is f l2.18p with r = 1. In 
this case, the ensuing nonlinear boundary condition also maintains the bulk symmetry u — )■ —u. 

Also note, the results are consistent with an interesting observation by Habibullin [8] who pointed 
out how in the classical sine-Gordon theory a general boundary condition could be obtained from 
a general Dirichlet condition (that is, one with uq ^ 0), using a combination of two carefully 
designed Backlund transformations. Explicitly, the combined effect of the two Backlund trans- 
formations is effectively identical to the Ghoshal-Zamolodchikov general boundary condition, 
namely, 

/n +\ ■ + ^o) , -1 • - '"o) 
Ux\}J,t) = Ksm. ^ h K sm ^ , 

where k is the parameter introduced by the Backlund transformations. It was pointed out also 
some time ago [13] that a type I defect is, from a classical point of view, a Backlund transfor- 
mation 'frozen' at the location of the defect. Therefore, the action of the defect in modifying 
a boundary condition is essentially the Habibullin construction and (12. 9p is the quantum field 
theoretic version of this. Bajnok and Simon [12] noted the same phenomenon, albeit stated 
differently, by adding a type I defect to a Dirichlet condition in the sinh-Gordon model. 

Using a type II defect to generate solutions of (12. 3p leads inevitably to solutions that are outside 
the Ghoshal-Zamolodchikov class since the explicit dependence of their coefficients on topological 
charge is unavoidable. One may wonder what boundary conditions might be responsible for 
these and one suggestion has been proposed with the boundary Lagrangian density (I2.20p . The 
most general solution (12. 8 p to the transmission Yang-Baxter equation in the form (12. 3p has 
been derived explicitly and a comprehensive explanation concerning the manner in which it is 
obtained using defects has been provided. One could imagine that repeatedly placing type I or 
type II defects near a Dirichlet boundary, or near a general Ghoshal-Zamolodchikov boundary, 
should generate all possible solutions to a (suitably generahsed) version of the compatibility 
relations (12. 3p containing additional pairs of topological charge labels (for example, see (12. 9p ). 
However, the types of classical boundary condition to which they might effectively correspond 
(if any) is unclear. One should bear in mind that any process of reduction (or 'fusing'), in 
which the additional sets of labels collapse to a single set, must yield the already known general 
solutions provided above, or a combination of them. On the other hand, the most general 
solution to (I2.10p (with a single pair of topological charge labels) is provided by (I2.12p (see 
also [28] for an alternative algebraic argument). Hence there is no other solution to (I2.10p that 
could be used in combination with a Dirichlet boundary condition to provide a more general 
solution to (12. 3p . It is worth recalling that boundaries with additional degrees of freedom have 
been considered previously for the sine-Gordon model by Baseilhac and Delius [29] and the 
corresponding reflection matrix found in [SU]- However,the latter represents a solution of a 
boundary Yang-Baxter equation that does not preserve the topological charge of the system. 
Hence it is not of the form (12. 3p . 

Finally, generalising to the complex affine Toda models, the analogue of the Dirichlet condition 
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is given in (11. 2p and it would be interesting to discover how the reflection matrices corresponding 
to these are related (if indeed they are) to the reflection matrices associated with the boundary 
conditions (ll.ip . The purpose of this note was to emphasise the role defects might play and 
these specific questions will be addressed in the future. 
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A Infinite dimensional solution to the reflection Yang-Baxter equa- 
tion 

Given the expression (12. 6 p for the reflection matrix the reflection Yang-Baxter equation becomes 
a collection of nonlinear recurrence relations for the coefficients. These are of four basic types, 
2-term, 4-term, 5-term and 6-term relations. It is also useful to define a^{xi,X2) = a(^2 ± 

^l), 6±(Xi,X2) =6(^2±^l). 

The 2-term relations involve just the off-diagonal elements and are: 

s±{a, xi)s±{a + 2, X2) = s±{a, X2)s±{a + 2, xi), 

s+(a, xi)s_(a -|- 2, X2) = s+(a, X2)s_(a -|- 2, xi), (A.l) 

implying s±{a,x) = f{x)s±{a), where f{x) depends on rapidity but not topological charge and 
the coefficients s±{a) do not depend on rapidity. 

There are a number of 4-term relations all of the following type, 

s+(a, X2) [a^b^ r+{a + 2, xi) — a~^b~ r+(a, zi)) 

= cs+(a, xi) (6"^ r+(a + 2, X2) - h~ r_{a + 2, X2)) , (A. 2) 

which simplify to 

/(X2) {a^h^ + 2, xi) — a^b" r^{a, xi)) 

= cf{xi) (6+ r+(a + 2, X2) - b~ r„(a + 2, X2)) , (A.3) 

and a number of 5-term relations of the form, 

s+(a, Xi) (a~'^a~ r^{a, X2) — b^b~ r^{a + 2, X2) — (?r^{a + 2, X2)) 

= cs+(a, X2) {a^ xi) — a~ r_(a + 2, xij) . (A. 4) 
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Again, using ( ]A.1|) . these simplify to, 

/(xi)((a"a"^ r+(a,X2) - r+(a + 2, X2) - r_(a + 2, 0:2)) 

= c/(x2) (a"*" r+(a;, Xi) — r„(a + 2, xi)) . (A. 5) 

Finally, there is a 6-term relation that reads 

c6~(r+(a,Xi)r+(a,X2) - r_(a, xi)r_(a, X2)) 

+ c6+(r+(a, xi)r_(a, X2) - r+(a, X2)r_(a, Xi)) 

+ a^6~(s+(a, xi)s_(a + 2, X2) — s+(a; — 2, xi)s_(a, X2)) = 0. (A. 6) 

The latter is the only nonlinear relation to supply a real constraint on the off-diagonal functions 
of a, s±{a). Making the ansatz /(x) ~ (x^ — 1/x^) and assuming r± are cubic in x and 1/x 
leads (using Maple) to the solution given in fl2.8p . In effect equations (1A.3P and (lA.Sp are solved 
first for r± and the last equation (1A.6P strongly constrains s±{a). 

More generally, it is useful to consider the zeros of /(x) and compare them with the zeros of 
r±{a,x). Suppose X2 is a root of /(x) and xi is not a root of /(x). Then, the right hand side 
of (1A.2P becomes 

b'^r^{a, X2) — b~r^{a, X2) = 0, 
and therefore, using the explicit forms of 6^, 

xj [x2r+(a,X2) + r_(a,X2)] - [x2 r_(a, X2) + r+(a, X2)] =0. 

Thus, since this must hold for any choice of xi, 

X2r+(a,X2) + r_(a,X2) = 0, X2r_(a,X2) + r+(a,X2) = 0, 

which together imply (1 — X2)r-|- = 0. Thus, either X2 is also a zero of r-|-(a,X2), or X2 = 1. 
In other words, the only zeros of /(x) not shared by the diagonal elements are located at the 
fourth roots of unity. 
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